跳到主要内容

Gram-Schmidt 正交化

阐述

  • q1q_1a1a_1 归一化得到
  • q2q_2a2q1q1a2a_2-q_1q_1^*a_2 归一化得到
  • q3q_3a3q1q1a3q2q2a3a_3-q_1q_1^*a_3-q_2q_2^*a_3 归一化得到
  • …… 

从这些就可以算出 RR 的系数。

算法:

这大约需要 2mn22mn^2 次浮点运算。但是,问题在于这样相当于不断用 Rj1R_j^{-1} 来右乘 AA,其中有可能 rjj0r_{jj}\approx 0,这样会导致 RR 是病态的。

所以,可以从另一个方面着手,用正交矩阵来左乘 AA,这样可以保证 κ(Qj)=1\kappa(Q_j)=1.

实例

性质

可以用于 QR 分解

相关内容

参考文献